Localization-delocalization transition of disordered d-wave superconductors in class CI 



O 

o 
o 

(N 

D 
tin 



Takahiro FukuitI 

Department of Mathematical Sciences, Ibaraki University, Mito 310-8512, Japan 

(January 24, 2000) 

A lattice model for disordered d-wave superconductors in class CI is reconsidered. Near the 
band-center, the lattice model can be described by Dirac fermions with several species, each of 
which yields WZW term for an effective action of the Goldstone mode. The WZW terms cancel out 
each other because of the four-fold symmetry of the model, which suggests that the quasiparticle 
states are localized. If the lattice model has, however, symmetry breaking terms which generate 
mass for any species of the Dirac fermions, remaining WZW term which avoids the cancellation can 
derive the system to a delocalized strong-coupling fixed point. 
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Dirty superconductors have attracted much interest, 
since they provide wider universality classes of disordered 
systems [Q-^. In particular, it is quite interesting to ask 
what the universality class of the disordered d-wave su- 
perconductors is 0-01 . A remarkable property of this 
unconventional superconductors is that near the band 
center, quasiparticle states can be described by Dirac 
fermions ||^. Such a description enables us, for exam- 
ple, to relate d-wave superconductors with the quantum 
Hall effect and to predict a new spin phase called spin 
quantum Hall fluid They should also produce the 
well-known effect of chiral anomaly ||l^,p| to c?-wave su- 
perconductors. Recently, Senthil et alf.] have studied 
disordered d-wave superconductors with spin rotational 
symmetry and reached the conclusion that all states are 
localized. One knows, however, that the WZW term due 
to chiral anomaly plays a crucial role in two dimensional 
critical phenomena Therefore, it is quite important 
to take the WZW term missing in [|| into account, or to 
answer the question why it vanishes if it does not exist. 

In this paper, we reconsider disorder effects on the 
quasiparticle properties of the d-wave superconductors in 
the class CI (those with spin rotational and time-reversal 
symmetries) using a replica technique. It is shown that 
each Dirac fermion associated with four nodes creates the 
WZW term. It turns out that they cancel each other and 
the resultant nonlinear sigma model suggests localization 
of the quasiparticles in dirty d-wave superconductors, as 
was shown by Senthil et al It should be stressed, 
however, that the WZW term is potentially realizable 
and the cancellation is accidental: It is due to the four- 
fold symmetry of the model. Therefore, if such symmetry 
is broken, the system flows to the strong-coupling fixed 
point described by the WZW model. 

Let us begin with a lattice Hamiltonian for singlet su- 
perconductors 



H 




A- 



(1) 



tji and 



We can choose real and symmetric matrices ti 
Aij = Aji taking account of the hermiticity as well as the 
spin-rotational and the time-reversal symmetries. In the 



absence of randomness, we choose the following parame- 
ters for a pure Hamiltonian with d-wave symmetry 



Ajj±£ - -Ajj±y - Ao, (2) 



where x — (1,0) and y — (0, 1). Moreover, consider in- 
troducing small terms Hi to break the CI symmetry, 

^j,j+x±y = Aj+x±y,j = ±^-^, Ajj lAi. (3) 

It will be shown momentarily that this parameter con- 
trols the localization-delocalization transition of the 
present model. 

The pure Hamiltonian Hq has four nodes, where gap- 
less quasi-particle excitations exist . Therefore, we can 
firstly take the continuum limit around the nodes of Hq 
and next incorporate the continuum expression of iJi, 
provided that Hi is small. The lattice operators are then 
described by the continuum slowly- varying fields near the 
band center as 

Cji/a ^ f''+^y'ipll{x) + r^'^'^^i^^^lix) 



+i-^-+^yiPll{x) + i^-^-^yip'^lix), 



2t/ 



(4) 



where a is a lattice constant, x — aj, and two kinds 
of lower indices of the field are referred to as spin, 
left-right (LR) movers, respectively, and upper index as 
node. Namely, the field variable lives in the space 
V ^C'^®C'^® C^. The pure Hamiltonian in the con- 
tinuum hmit is then Ha = J (Pxij}'^ {Ho+'Hi)%l^ with 



(x ^ y) 



(5) 



where the coordinates have been transformed as 



x,y 



. The explicit matrix in Eq. (||) denotes the node 
space, and matrices 7^ belong to the other space of V , 
calculated initially as 71 = vf'^2®<^5 and 72 = WAl2®cri, 
where vp = 2^/2toa and wa = 2V2Aoa. It may be more 
convenient to choose 



1 



7l = VFI2 ® 0-2, 72 



'VAh 



(6) 



via suitable rotation in LR-space of V. In this basis, Hi 
is given by 



Hi = 







ml2 (8) 0-3 



(7) 



where m = Ai. Namely, -ffi yields asymmetric mass 
term in the continuum Hamiltonian. For the time be- 
ing we neglect it, but it will be shown that vanishing m 
leads to localization whereas finite m drives the system 
to delocalization. 

The spin-rotational and time-reversal symmetries of 
the lattice model translate, respectively, into the con- 
tinuum model as 



n 
n 



c - 
T 



1(72 ® O"! 

= I2 o^s ( 



)l2, 

12, 



(8) 



where t means the transpose. The total Hamiltonian 
density is given hy H = Ho + Hd, where Tid is disorder 
potential satisfying Eq. (^). It is stressed that we take 
account of all kinds of disorder potentials satisfying Eq. 
(^. This implies that the summation with respect to i, j 
in Eq. should be over on-site, nearest-neighbor, and 
diagonal-second- neighbor pairs. Namely, we can achieve 
"maximum information entropy" for the Dirac Hamil- 
tonian with the symmetries (|^) when we introduce not 
only on-site disorder potentials but also disordered hop- 
ping and off-diagonal pairing for the lattice model in Eq. 
(|l|). Actually Ludwig et al have derived generic Dirac 
Hamiltonian for the integer quantum Hall transition con- 
sidering lattice model in a similar situation It may 
be straightforward to explicitly calculate disorder poten- 
tials but tedious to average over them, because there are 
no less than twenty independent potentials satisfying Eq. 

. The key point for the ensemble average is that if the 
model has "maximum entropy" , we can use the technique 
developed by Zirnbauer Q. 

To study one-quasiparticle properties of the model, we 
introduce the Green function Gaa'{x, x'; ie) — {x, a\{ie ~ 
H)~^|a;', a'), where index a denotes the set of spin, LR, 
and node species in the space V and |x,a) = ■0l(a^)|O). 
Especially, we need 



a 

K{x,x') = Gaa>{x, x';ie)Ga'a{x', x; -ie) (9) 

a, a' 

to compute the DOS and the conductance of the quasi- 
particle transport jl^ . Some notations are convenient 
for this purpose. The introduction of replica for the 
field ip enables us to express the generating functional of 
these Green functions by path integrals: ipa — > V'aa and 
ijj* ipaa^ where a and a are indices denoting V and the 
replica space Wi- = C", respectively. The fields ijj and 



tp* have been converted into matrix fields, which makes 
it simpler to define an order parameter field. It should 
be stressed that the fields -0 and ip* are completely inde- 
pendent variables. Lagrangian density is then described 
symbolically as £ = — trw,,?/'^ (*e~H)V'j where tr^^ is 
the trace in the replica space and the summation over 
the indices a of the V space is implied according to the 
rule of the matrix product. It should be noted again that 
'0^ is independent of ip. Moreover, we introduce an aux- 
iliary space Q to reflect the symmetries in the V space 
(0) to an auxiliary field introduced later [See Eq. (fsl) 



Wr ^ W = W,(g) Wa with VFa = «) which are 
associated with the spin-rotational and the time-reversal 
symmetries, respectively. Fermi fields are now denoted 
by and ^'^q. One of simpler choices is 



i+ \ ,7;, - ,7,u 
0- 



(10) 



where 7± — (1 ± T)/2 serves as a projection operator ( 
T+-fT_ = I2«)l2«)l2, Tl^T±, and T+T_ = ) into each 
chiral component of Sp(n)xSp(n) symmetry, as we shall 
see later, and — \ (01, — iV'2) with 0i,2 = ± iC^^i)* . 
The newly introduced fields ^ and ^ are subject to |^, 



Matrices 7 and tt are defined in the W space by 



(11) 



7 

TT 



' icT2 ® I2 
) I2 ® era. 



70 ® I2, 



(12) 



The identity trH/(«ew^'«' - ^H"^) = tYw,4>\ie - n)i^, 
where uj = 1„ ® 0-2 cri, leads to the generating func- 



tional Z ^ J V^V^e 



with 



S 



(13) 



Assume that disorder potential Tid obeys the Gaus- 
sian distribution P[Hd] — exp^— ^trvH^^ Then en- 
semble average over disorder is quite simple. The pro- 
cedure is as follows Q|: The disorder potentials are inte- 
grated out by using the identity —j^ivv'H^+tvv'Hd'^'^ = 
-^(tryHd - g^**)^ + ftry(1'$)2. It turns out that 



the integration over Tid is automatic because vj/xp sat- 
isfy the same symmetries as those of Tid due to Eq. 
([III). This is actually a merit to consider the disorder 
potentials with maximum entropy. If disordered hop- 
ping and off-diagonal pairing of the lattice model are ne- 
glected and on-site disorder potentials are merely taken 
into account, some other conditions should be imposed 
on Tid and therefore on the fields 4" and 5*. Now we 



2 



have interaction terms of fermions due to ensemble aver- 
age. Note the identity try(^'$)2 = -trwi^^)'^- Then, 
the four fermi interactions are decoupled via auxiliary 
matrix (order parameter) field defined in the W space. 
To be concrete, add the following term into the action, 
j^tiwiQ + g'^'^ — i^)"^, which is actually a constant af- 
ter integration over Q. Then we reach an effective La- 
grangian density. 



M 



g-Mg\, i^±^M±, ^±^g±i^±- (20) 



— (Q2 „ 2ieiuQ) + Q** - *Ho* 



(14) 



Here we have set J = for simplicity. Notice that the 
anti-Hcrmitian auxiliary field Q = —Q^ is subject to 



Q = -lQ'l~\ Q = -7rQTT-\ 
The solution of these equations is 

-Aft 

M 



(15) 



(16) 



with a condition M = 7oM*7g ^, where the explicit ma- 
trix in the above equation denotes the time-reversal space 
of M^. 

The Lagrangian (|lj) has G=Sp(n)xSp(n) symmetry. 
To see this, let us consider the transformation Q — > 
gQg~^ which keeps the symmetry relations ([l5|). It turns 
out that g should satisfy 7 — g^g^ and rgr^^ = 5 as well 
as gg^ = 1, and therefore g is explicitly given by 



g 



g+ 



(17) 



where g± £ Spin) is a 2n x 2n matrix in the replica 
and the spin space of 14^. So far we have derived the 
Lagrangian as well as its symmetry group. To inte- 
grate out the fermi fields, it may be convenient to use 
the notations where the fermi fields ^' and "if are col- 
umn and row vector, respectively, as usual. Using Eq. 
(|l5|), fermion part of the Lagrangian (|lj) is rewritten as 
Cf = $(7^o«)l-l(g)g*)4^= *l(8)7(Ho®l+l(8)(3)l«)7t«'. 
Transform the fields as ^ (g) 7'l' and * ^ 1 (Ki 7*. 
Then the Lagrangian is given by, in terms of the fields 
M and ip±, 

£= itriv^ [M^M - e7o(A/ - AP)] + Cfi + Cf2, (18) 

g 

where W-^s is the replica and the spin space of W and C-pj 
describes the Lagrangian of the jth-node fermion defined 
by 

£fi = 4j\i^ijX + - tPIM^'tPI + rplMiPl, (19) 

and = >Cfi(1 — > 2,a; y). Here and hereafter, the 
identity matrices such as those in/) ^ 1 and 1 (g) Af are 
suppressed. The transformation laws of M and ^± fields 
are 



We have used a bit complicated basis for W in the defi- 



nition of and ^ in Eq. ( 10|) , since Q and g become sim- 
pler in this basis. On the other hand, it may cause a dif- 
ficulty in computing the saddle points. We have known, 
however, from the e-term in Eqs. (|l3|) and (^4|) that oj 
serves as a "metric" in the extended auxiliary space. Usu- 
ally we may choose a basis of the space W with a diagonal 
metric, owing to which we can assume that Q is also diag- 
onal on the saddle points. In the present case, therefore, 
it is natural to assume that Q should have the same struc- 
ture as 07, and hence Mq = vjq with real diagonal matrix 
V — diag(ui, • • • , w„). Then the variation with respect 
to V after the integration over fermi fields tells that the 
saddle points are given by Va = vq ^ Aexp{~Trvpv^/ g), 
where A is a ultraviolet cut-off. This solution gives rise 
to an exponentially small density of state at the band- 
center. Now it is easy to identify the saddle point mani- 
fold as H=Sp(n): The chiral transformation g = {g+, g-) 
in Eqs. ( |l7|) and ( po|) is divided into two types. One is 
5v — (ffiiSi) under which Mq is still invariant, and the 
other is ga — (52, si) under which Mq is no longer invari- 
ant. 

Nonlinear sigma model is derived as small fluctua- 
tion around the saddle point manifold H by consider- 
ing type local Sp(n) rotation. Let us parameterize 
M = CH^ = HU with C e Sp(n), H = CH^\ and 
U = C^. The field ^ describes the massless fluctua- 
tion around the saddle point manifold H, whereas H de- 
scribes massive longitudinal modes. In what follows, we 
take only the leading order for the latter mode, setting 
H = A/o(= H). It should be noted that nonlinear sigma 
model on G/H has a global G symmetry as well as a hid- 
den local H symmetry. Though the field ^ itself is not 
invariant under local H transformation, the composite 
field U is invariant. 

To derive an effective action for the transverse mode, 
let us come back to the Lagrangian (|l9|), since we should 
be careful in the integration over fermi fields. It is practi- 
cal to firstly integrate out the fermi field of node-1. Then, 
the contribution from node-2 can be obtained by replac- 
ing x <-> y. To carry out the former integration, make 
the transformation x+ — '4'+U'^ and x+ = whereas 
X- = V'-ln ® f2 and X- = In ^ cr2ip-- The Lagrangian 
( p9| ) is then converted into 



£fi - x\iPx\- + X-i^X- + ^vo {x\x- + X-X\-) 



(21) 



-L^ with L^^Ud^^W. 



where D^j_ is defined by = 9^ 
It is convenient to scale x = vpx' and y = v^y' working 
in the node-1 sector. Integration over the fermi field of 
node-1 yields 



IVQ 



i{-idi - 82) 



i{iDi - D2 



IVQ 



(22) 



3 



where the identity matrix I2 belongs to the spin space 
of V, Vi means the node-1 sector of V, the derivatives 
are with respect to the scaled coordinates x' and y', 
and ^i{U) is the Jacobian for the chiral transformation, 
which can be calculated by using the Fujikawa method 
| p3[ . We simply present the final answer Zpi ~ e^'^\ 
where the effective action Si associated with the nodc-1 
fermion is composed of the principal chiral action of U 
with a coupling constant X = 4tt and of the WZW term 

rwzw = -^ [ d'^xe^^rti-w^d^UU^dMU^drUU^. (23) 

IZTT J 

Therefore, integration over fermion of node-1 actually 
yields the WZW term. 

Next let us compute the contribution from the node-2. 
The procedure is the rescaling x' — x/vp and y' = yjvi^ 
and the exchange x ^ y. It should be noted that the 
WZW term is invariant under the rescaling but is odd 
under the exchange, and hence it cancels out. The total 
effective action ends up with 



d xirw . 



kV 



WZW 



(24) 



with 4 — and fc = 0. This is just the action 

A 47r vfva •' 

derived by Senthil et al. Although the WZW term dis- 
appears in the ordinary c?-wave superconductors, we are 
tempted to expose it, since the WZW term exists po- 
tentially. This is indeed possible: If the lattice model 
includes the symmetry breaking term (^ and therefore 
the Dirac fermion for the node-2 has a mass (0), we can 
neglect the node-2 fermion in the lower energy than the 
mass gap. In this case, we have the same action (^) but 
with T = -r- and fc = 1 in the scaled coordinates x' and 

y'- _ 

The renormalization group equations of the action (|24| ) 
are calculated at the one-loop order as 



d In L 47r 

rfe _ A 
dlnL V Svr 



47r 



(25) 



where d = 2, e — d — 2, and the replica limit n ^ has 
been taken. In the case where m = 0, we reach the same 
conclusion as Senthil et al. Namely, one-quasiparticle 
states are localized, since the spin conductance against 
weak magnetic fields is related with A as ct — 2/ (ttA) , 
which is calculated Jo via diffusion constant in the dif- 
fusion propagator (P). It should be stressed that the 
cancellation of the WZW term is due to the four-fold 
symmetry of the d-wave Hamiltonian. Actually, the la- 
tent WZW term can emerge via symmetry breaking mass 
term for the Dirac fermions, and in that case the cou- 
pling constant A flows to the strong-coupling fixed point 



value Ac = 47r/fc. This fixed-point is conformal invariant 
described by Sp(n) WZW model. From the scaling di- 
mension of the energy in Eq. ( |25| ) , it turns out that the 
density of state near this fixed-point obeys the scaling 
law 



(26) 



Therefore, if the pure model has the breaking term (|^), 
we suggest that p{E) = Ei . 

In this paper, we have taken the breaking term of type 
(^) into account. More detailed phase diagram will be 
published elsewhere. Numerical check of the present con- 
jecture does not seem difficult. It is quite interesting to 
expose the hidden WZW term which exists potentially 
but conceals itself in the four-fold symmetry of the d- 
wave superconductors. 
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